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ABSTRACT 

Recent spectroscopic observations in the Milky Way suggest that the chemically defined thick disk 
(stars with high [a/Fe] ratios and thus old) has a significantly smaller scale-length than the thin disk. 
This is in apparent contradiction with observations of external edge-on galaxies, where the thin and 
thick components have comparable scale-lengths. Moreover, while observed disks do not flare (scale- 
height does not increase with radius), numerical simulations suggest that disk flaring is unavoidable, 
resulting from both environmental effects and secular evolution. Here we address these problems 
by studying two different suites of simulated galactic disks formed in the cosmological context. We 
show that the scale-heights of coeval populations always increase with radius. However, the total 
population can be decomposed morphologically into thin and thick disks, which do not flare. We 
relate this to the disk inside-out formation, where younger populations have increasingly larger scale- 
lengths and flare at progressively larger radii. In this new picture, thick disks are composed of the 
imbedded flares of mono-age stellar populations. Assuming that disks form inside out, we predict 
that morphologically defined thick disks must show a decrease in age (or [a/Fe] ratios) with radius 
and that coeval populations should always flare. This also explains the observed inversion in the 
metallicity and [a/Fe] gradients for stars away from the disk midplane in the Milky Way. The results 
of this work are directly linked to, and can be seen as evidence of, inside-out disk growth. 

Subject headings: Galaxy: evolution — Galaxy: abundances — galaxies: kinematics and dynamics 


1. INTRODUCTION 

The formation of galactic thick disks has been an im- 
por tant topic ever since their disco very in external galax¬ 
ies (lBursteinlll979t lTsikoudilll979D and in the Milky Way 


([Gilmore and Reidf [T983). A number of different mecha¬ 
nisms have been pro posed fo r th e format ion of thick disks 
(see Introduction in [Minchev et al.1 12012| and references 
therein), but it is generally believed that they are the 
oldest disk components. 

Stellar disk density decomposition into thinner and 
thicker components in external edge-on galaxies have 
found that the thin and thic k disk components h ave 
comparable scale-lengths (e.g.. lYo achim a nd Da lcanton 
120061: IPohlen et ahll2007t IComeron et al.l 12012 1. While 
this is consistent with results for the Milky Way when 
similar morphologically ( or structural) definition for the 
thick disk is use d (e.g., iRobin et al.l 119961 : lOihal [20071 
[Juric et ahl l2008i l. it is in contradiction with the sig¬ 
nificantly more centrally concentra ted o l der or [q/Fe]- 
enhanced stellar populations (e . g.. iBensbv et all 120111 : 
ICheng et al.l l2012at iBovv et~ahl I2012D . This apparent 
discrepancy may be related to the different definition of 
thick disks - morphological decomposition or separation 
in chemistry. 

Moreover, w h ile ob ser ved disk s do not flare 
(Ivan der Kruit fc Searlel Il982t Ide Griisl 119981 : 
IComeron et all 1201 lh . numerical simulations sug¬ 
gest that flaring cannot be avoided due to a range 
of different dynamical effects. The la rgest source is 
probably satellite-disk interac t ions l'e.g..lBournaud et al.1 
l2009t iKazantzidis et al.l 120081 : IVillalobos fe Helmil 120081) 


1 Leibniz-Institut fur Astrophysik Potsdam (AIP), An der Stern- 
warte 16, D-14482, Potsdam, Germany 

2 Max-Planck-Institut fur Astronomie, Konigstuhl 17, D-69117 
Heidelberg, Germany 


shown to increase an initially constant scale-height by 
up to a factor of ~ 10 in 3-4 disk scale-leng ths. S t udying 
preassembled N-body disks, I Minchev et al.l ( i2012n 
showed that purely secular evolution (in the absence 
of external perturbations) also causes flared disks, due 
to the redistribution of disk angular momentum. This 
was linked to the accumulation of outward/inward 
migrators in the outer/inner disk (with vertical actions 
larger/smaller than those of non-migrators), which 
caused vertical thickening/contraction and thus flar¬ 
ing. In addition to merger perturbat ions a nd ra dial 
migration, misaligned gas infall (IScannapieco et all 
2009; Roskar ct al. 2010) and reorien tation of the disk 
rotation axis (lAumer & Whitel 12013 ! can also produce 
disk flaring. 

In this Letter we show that the above outlined contra¬ 
dictions can be reconciled when realistic galactic disks 
growing in a cosmological environment are considered. 


2. SIMULATIONS 

We study two simulations following the formation of 
disks in the cosmological context, which use two distinct 
simulation techniques. 

The first model, Modell, is a gas dynamical simula¬ 
tion in a cosmological c ontext using a sticky parti cle al¬ 
gorithm, as described in IMartig et all (|2009L [2?)T2h . with 
150 pc spatial and 10 4-5 M 0 mass resolution. This is a 
barred late-type galaxy with stellar mass 4.3 x 10 1o Mq 
and a disk scale-length hd = 5 kpc, estimated from 
all stars in the range 4 < r < 12 kpc, \z\ < 2 kpc. 
Further deta il s abo ut this simulation c an be found in 
IMartig et al.l (1201411 (model gl06) and iMinchev et al.1 
(120131b 

The second model, Model2, is a full cosmological zoom- 
in hydro simulation, using initial conditions from one 
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Fig. 1.— Scale-height fits for Modell. Top: The vertical density 
profiles of mono-age populations (a typical one shown here) can 
be fitted well by single exponentials at all radii. Different colors 
correspond to different radial bins, as indicated by the color-bar. 
Distances are in units of the disk scale-length, hd . Bottom: When 
the total stellar population is considered, for all radii a sum of 
two exponentials is required for a proper fit. Model2 gives similar 
results. 


of the Aquarius Projec t haloes (iSprineel et al.l 120081 : 
IScannapieco et all i2009T) . The technique used here is 
Tree-PM SPH with 300 pc spatial and 4.4 x 10 5 M Q 
mass resolution. Model2 is a non-barred galaxy with 
stellar mass 5.5 x 10 1 o Mq and a disk scale-length hd = 
4 kpc, estimated in the same spatial range as for Modell. 
Furt her deta ils ab out this simulation can be found in 
lAumer et al.l (|2013l h their model Aq-D-5. 

Both simulations form initial central components dur¬ 
ing an early epoch of violent merger activity. Gas-rich 
mergers supply the initial reservoir of gas at high red- 
shift and merger activity decreases with redshift, sim¬ 
ilarly to what is expected for the Milky Way. This 
inside-out disk formation results in centrally concen¬ 
trated older stellar pop ulations for both si mulations (see 
IMartig et al.l 120141 and lAumer et al.l I2014D . The general 
formation and evolutionary behavior of our disks is sim¬ 
ilar to ma ny recent simulations in the cosmologi cal con¬ 
text fe.g.. lBrook et al.ll2012l : IStinson et al.ll2013h . 

We work in units of disk scale-lengths, hd , which are 
estimated using the total stellar population at the final 
simulation time. 


3. RESULTS 


3.1. To flare or not to flare 

Similarly to IMartig et al.1 (120141 ) , we decomposed the 
stellar disks into mono-age populations, i.e., narrow bins 
of age, where we used Aage = 1 Gyr. It was found that 
single exponentials provided good fits to the column den¬ 
sity in the vertical direction for all age b ins and bot h 
Mode ll and Model2, in agreement with IMartig et al.1 
(12014T) . Assuming that chemistry can be a proxy for age, 
this is also consistent with stud ies o f mon o-ab undance 
populations in the Milky Way dBovv et al.ll2012l ). The 
top panel of Fig. [T] shows single exponential fits for a typ¬ 
ical mono-age population of Modell, in eight radial bins 
in the range 1 < r < 5 disk scale-lengths, hd- Radial 
bins of width 0 .6hd were used. It can be seen that the 
fits are very good for all radii. In contrast, the bottom 
panel of Fig. [T| shows that a sum of two exponentials is 
required for a proper fit when considering the total stellar 
population. 

It is clear already from the top panel of Fig.[T]that sig¬ 
nificant flaring exists for the particular mono-age popula¬ 
tion shown (i.e., the vertical density profile flattens with 
increasing radius). However, only milder variations with 
radius are seen for the total population, when decom¬ 
posed into thin and thick disks (Fig. [TJ bottom panel). 

In the top row of Fig. [2] we plot the scale-height vari¬ 
ation with galactocentric radius, r, in the region 1-5 
disk scale-lengths, hd, for both Modell (left panel) and 
Model2 (right panel). Forty overlapping mono-age popu¬ 
lations of width 1 Gyr are shown, covering the age range 
in each disk, as indicated in the color bar on the right of 
each panel. Both the radius and scale-height, h z , are in 
units of hd- It can be seen that, indeed, for both models 
significant flaring is present, which increases for older co¬ 
eval populations. For Modell the increase in scale-height 
with radius for stars of age ~ 8 Gyr is about a factor of 
four in four scale-lengths. The flaring decreases by about 
a factor of two for stars with age < 2 Gyr. Even more 
significant flaring is seen for Model2. 

In contrast to the flaring found for all mono-age popu¬ 
lations, the thin and thick disk decomposition of the total 
stellar population including all ages, results in no appar¬ 
ent flaring. This is shown by the triangle and square 
symbols overlaid in the top row of Fig. [2] 

We discuss the interpretation and implications of this 
interesting result in the following sections. 

3.2. Interpretation 

What is the reason for the flaring of mono-age disks? 
In numerical simulations flaring is expected to result 
from a number of mechanisms related to galactic evo¬ 
lutionin a cosmological context (e.g., [B ournau d et al.1 
20091: iKazantzidis et al.l 120081 : lYillalobos &; Helmil 120081 : 
Aumer fe Whitell2013li . Even in the absence of environ¬ 
mental effects, flaring is unavoidable due to secular evo¬ 
lution alone (radial mi grati on caused by spirals and/or a 
central bar, iMinchev et al.ll20l2lh It should be stressed 
here that, while migration flares disks in the lack of ex¬ 
ternal perturbations, durin g satellite-disk interact ions it 
works against disk flaring (Min chev et al.ll2014bli . Yet, 
this is not sufficient to completely suppress the flaring 
induced by orbiting satellites, as evident from the top 
row of Fig. [2] This suggests that external effects are 
much more important for the disk flaring we see in the 
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Fig. 2.— Top: Variation of disk scale-height, h z , with galactic radius for Modell (left) and Model2 (right). Color lines show mono-age 
populations, as indicated. Overlapping bins of width Aage = 1 Gyr are used. Overlaid also are the thin (triangles) and thick (squares) 
disks obtained by fitting a sum of two exponentials to stars of all ages (as in Fig. [T] bottom). No significant flaring is found for the thin 
and thick disks. Middle: Disk surface density radial profiles of mono-age populations. Older disks are more centrally concentrated, which 
explains why flaring diminishes in the total population. Also shown are the surface density profiles of stars close to (triangles) and high 
above (squares) the disk midplane. The thickened disk component extends farther out than the thin one. Bottom: Variation of mean age 
with radius for samples at different distance from the disk midplane, as indicated by the color bars. Slices in |^| have thickness A|z| = 0.1 hd- 
Overlaid are also the age radial profiles of stars close to (triangles) and high above (squares) the disk midplane. Age gradients are predicted 
for both the (morphologically defined) thin and thick disks. 


simulations. Because the mass and intensity of orbiting 
satellites generally decreases with decreasing redshift, so 
does the flaring induced. It can be expected that at a 
certain time secular evolution takes over the effect of ex¬ 
ternal perturbations. 3 


Resolution and sub-grid models for gas physics may set 
an unrealistically high lower lim it for the vertical v elocity 
dispersion of young stars (e.g., iHouse et al.ll201lh . The 
fact that the youngest populations in our simulations do 
not flare (dark-blue curves in top panels of Fig. [2| seems 


3 IMinchev et alj Il2014af ) suggested that this time could also be of narrow metallicity samples, 
inferred from the shape of the [a/Fe]-velocity dispersion relation 
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to indicate that nu merical issues are not important for 
our result (see also lMartig et al.ll20l4 Sec. 5). 

What is the reason for the lack of flaring in the to¬ 
tal disk population? In an inside-out formation scenario, 
the outer disk edge (where flaring is induced) moves pro¬ 
gressively from smaller to larger radii, because of the 
continuous formation of new stars in disk subpopula¬ 
tions of increasing scale-length. At the same time the 
frequency and masses of perturbing satellites decreases. 
Because of the inside-out disk growth, which results in 
more centrally concentrated older samples (see Fig. [2l 
middle row), the younger the stellar population, the fur¬ 
ther out it dominates in terms of stellar mass. Therefore, 
the thickened disk component results from the imbedded 
flares of different coeval populations, as seen in the top 
row of Fig. [2] 

From the above discussion, it can easily be seen that 
inside-out disk growth would always have the effect of 
suppressing disk flaring in the total stellar population, 
if not completely removing it. We produced the equiv¬ 
alents of Fig. [2] for other disks from the two simulation 
suites considered here. The flaring was always strongly 
reduced when the total population was considered, sug¬ 
gesting that our results are generic. 


for age fe.g.. lMatteucci fe Brocatoltl99dlChiappini et al.l 

Horn 

Negative [a/Fe] g radie nts at large \z\ were r eported by 
IBoeche et al.l (1201311 and iBoeche et al.1 (12014D in dwarfs 
and gia nts, respectively, stud yi ng data from the R AVE 
survey dSteinmetz et alll2006f) . lAnders et al.l (120141) and 
lHavden et al.l (|2014ll also showed a transition from a 
weakly positive (or flat) to a negative radial gradient 
in APOGEE [a/M] measurements, as the distance from 
the Milky Way disk midplane was increased. Similar 
inversion of [a/Fe] gradients ha s been also found in 
recent chemo-dyna mical models dMinchev et al.l l2014bl : 
iRahimi et al.112014 : Miranda et al. 2015, submitted). 

Flaring of mono-age disks can result also in the 
inversion of metallicity gradients with increasing dis¬ 
tance from the disk midplane, as younger metal-rich 
stars in the outer disk can reach high vertical dis¬ 
tances. Observationally, this inversion has been found 


SEGUE - Chens et all [2012b 

RAVE - Boeche et al.1 

20131 APOGEE - Anders et al 

2014) and in simulations 

(Minchev et al.l 2014a; Rahimi et al.l 20141). 


All of the above results are in agreement with the de¬ 
cline in age with radius found in this work. 


3.3. Radial age gradients at high distance from the disk 

midplane 

From the top panels of Fig. [2] it can be expected that 
morphologically defined thick disks, as in observations 
of edge-on galaxies, will show a decrease in age with in¬ 
creasing disk radius. To quantify this, in the bottom row 
of Fig. [2] we show the mean age variation with radius for 
stellar samples at different distance from the disk mid¬ 
plane, as indicated by the color bar on the right. Slices 
in \z\ have thickness A\z\ = 0.1 hd- 

We also plot the mean age of samples close to (trian¬ 
gles) and high above (squares) the disk midplane, to be 
associated with thin and thick disks defined morpholog¬ 
ically. The two vertical ranges considered, \z\/hd < 0.1 
and 0.3 < \z\/hd < 1.0, for the Milky Way would corre¬ 
spond to \z\ < 0.3 kpc and 0.9 < \z\ < 3.0 kpc, using a 
disk scale-length hd = 3.0 kpc. 

A decrease of age with radius is found for the above 
morphologically defined thin disk in both models. How¬ 
ever, a significant negative age gradient is also seen in 
the morphological thick disks. Over the radial extent of 
4 scale-lengths shown in Fig. [2] we find a drop in mean 
age of ~ 3 Gyr and ~ 4.5 Gyr for the thick disks of 
Modell and Model2, respectively. 

4. IMPLICATIONS FOR THE MILKY WAY 

We can already look for evidence of our prediction in 
Milky Way observations. 

4.1. Inversion in radial [a/Fe] and [Fe/H] gradients 
away from the disk midplane 

As discussed above, our results suggest that even at 
distances high above the disk midplane, mean stellar age 
should decrease with galactic radius. As good age esti¬ 
mates for large stellar samples are currently unavailable, 
we can instead look at measurements of stellar [a/Fe] ra¬ 
tios. Because stars of different masses release chemical el¬ 
ements to the interstellar medium on different timescales, 
abundance ratios, such as [a/Fe], can be good proxies 


4.2. Flaring of younger stellar populations 

Studying APOGEE data, dNidever et al.l[20T4h showed 
that, at 1 < |z| < 2 kpc from the disk plane, the high- 
[a/Fe] sequence dominates the inner disk (5 < r < 
1 kpc), while in the outer disk (9 < r < 11 kpc) only 
the low-[a/Fe] sequence is present (see their Fig. 10). 

Flaring of younger stellar populations in the Milky 
Way have been reported in severa l works (e.g., 
(iFeast et al.l 12014 iKalberla et al.1 12014 iCarraro et al.1 
I2015D . 

All these results are in agreement with the idea of 
mono-age flaring populations suggested here, which al¬ 
low chemically- or age-defined thin-disk stars to reach 
large vertical distances in the outer disk. 


4.3. Extended morphologically defined thick disks 

We are suggesting here that inside-out disk for¬ 
mation inevitably gives rise to centrally concentrated 
older/[a/Fe]-enhanced populations, but at the same time 
- to comparable scale-lengths for the morphologically de¬ 
fined thin and thick disks. It can be seen in the middle 
row of Fig. [2] that this expectation is borne out for both 
our models. While inside 3-4hd mono-age populations 
become steeper with increasing age, the surface density 
profiles (of stars of all ages) close to (triangles) and high 
above (squares) the disk midplane run mostly parallel. 
Moreover, beyond 3-4 hd the thick disk profiles become 
flatter than the thin ones, i.e., extend further out. 

A morphologically defined thick disk in the Milky Way 
with a scale-length l arger than that of the thin disk has 
been reported by iRobin el al . 199(il using a compilation 
of d ifferent data sets, lOjlral (l200lD using 2MASS data, 
and l.luric et all (l2008h using SDSS data. This is well 
in line with external observations and explained by the 
idea presented here, resolving the apparent disagreement 
with the centrally concentrated [a/Fel- enhanced popula¬ 
tions in the Galax y reported bv lBensbv et al.l (l2011h and 
iBovv et all (|2012D . 
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5. DISCUSSION 

We showed that in galactic disks formed inside-out, 
mono-age populations are well fitted by single exponen¬ 
tials and always flare. In contrast, when the total stellar 
density is considered, a sum of two exponentials is re¬ 
quired for a good fit, resulting in thin and thick disks, 
which do not flare. This can be explained by the in¬ 
crease in scale-length of younger populations flaring at 
progressively larger radii. This scenario resolves the ap¬ 
parent contradictions that chemically (or age) defined 
thick disks are centrally concentrated, but structurally 
thick populations in both observations of edge-on galax¬ 
ies and in the Milky Way, extend as much as, or even 
beyond, the thinner component. 

The result that flared mono-age populations give rise 
to non-flared thin and thick disks is directly linked to 
the disk inside-out formation. If galactic disks formed 
with a constant scale-length as a function of time, then 
the flares of all mono-age groups would simply add up, 
making significant flaring of the total stellar population 
unavoidable. This is not seen in observations of edge-on 
galaxies. Our results, thus, land support to this widely 
accepted idea of inside-out formation of galactic disks. 

We found similar results in two simulation suites, rep¬ 
resented here by Modell and Model2, using radically dif¬ 
ferent simulation techniques. This suggests that our re¬ 
sults are generic and should be found in all cosmological 
simulations, where inside-out disk growth takes place. 

While flaring has been seen as a general characteristic 
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of disk heating by mergers, this conclusion can be linked 
to the initial conditions of preassembled disks typically 
used in such studies (allowing for a systematic explo¬ 
ration of parameter space). However, the results of this 
work suggest that heating by mergers can produce thick 
disks which do not flare, by considering the more realis¬ 
tic scenario of inside-out growing disks in a cosmological 
environment. It also implies that flaring in the total stel¬ 
lar population in observations of external galaxies cannot 
be used as a measure of the amount of mergers suffered 
by the host disk. Instead, flaring in narrow age groups 
should be searched for. 

The phenomenon described here results from the in¬ 
terplay among the intensity of inside-out formation and 
the masses of infalling satellites and their frequency as 
functions of redshift. Further work is necessary to find 
observables which can differentiate between these pro¬ 
cesses. Our results may then form the basis for future 
high-precision surveys (e.g., Gaia in the Milky Way) to 
constrain disk formation mechanisms by quantifying the 
amount of mergers and gas inflow as a function of cosmic 
time. 
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